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Abstract 

In economics, insurance and finance, value at risk (VaR) is a widely used 
measure of the risk of loss on a specific portfolio of financial assets. For a 
given portfolio, time horizon, and probability a , the 100a% VaR is defined 
as a threshold loss value, such that the probability that the loss on the port¬ 
folio over the given time horizon exceeds this value is a. That is to say, it 
is a quantile of the distribution of the losses, which has both good analytic 
properties and easy interpretation as a risk measure. However, its extension 
to the multivariate framework is not unique because a unique definition of 
multivariate quantile does not exist. In the current literature, the multivari¬ 
ate quantiles are related to a specific partial order considered in R”, or to a 
property of the univariate quantile that is desirable to be extended to R". In 
this work, we introduce a multivariate value at risk as a vector-valued direc¬ 
tional risk measure, based on a directional multivariate quantile, which has 
recently been introduced in the literature. The directional approach allows 
the manager to consider external information or risk preferences in her/his 
analysis. We have derived some properties of the risk measure and we have 
compared the univariate VaR over the marginals with the components of the 
directional multivariate VaR. We have also analyzed the relationship between 
some families of copulas, for which it is possible to obtain closed forms of 
the multivariate VaR that we propose. Finally, comparisons with other al¬ 
ternative multivariate VaR given in the literature, are provided in terms of 
robustness. 


1 Introduction 

Value at risk (VaR) has become a benchmark for risk management which is defined 
as the threshold quantity that does not exceed a certain probability level which is 
considered to be dangerous. It is commonly implemented by investment banks to 
measure the market risk of their asset portfolios. Although (VaR) has been broadly 
criticized from the work of [Artzner et al. (1999)] since it does not verify the di¬ 
versification property, it has also been defended by [Heyde et al. (2009)] for its 
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robustness. For univariate risks, the VaR is simply the ct—quantile of the loss dis¬ 
tribution function. Thus, the VaR is a risk measure easily interpretable, and it still 
remains the most popular measure used by risk managers. Unfortunately, a unique 
definition of multivariate VaR is more complicated because there are different pos¬ 
sible definitions of multidimensional quantiles that try to generalize some desir¬ 
able properties of the univariate quantile. For instance, the proposals given by 
[Koltchinskii (1997)] of multivariate quantiles as inversions of mappings, multi¬ 
variate quantiles in terms based on norm minimization as in [Chaudhuri (1996)], 
multivariate quantiles as level-sets given by [Fernandez-Ponce and Suarez-Llorens (2002)], 
multivariate quantiles based on depth functions developed in [Serfling (2002)], and 
finally, multivariate quantiles based on projections as in [Fraiman and Pateiro-Lopez (2012)], 

[Hallin et al. (2010)], [Kong and Mizera (2012)]. 

Currently business and financial activities generate data for which it has been 
shown that it is insufficient to consider single real-value measures over marginal 
aspects, in order to quantify risks jointly associated to the data. For instance, one 
of the drawbacks detected in the global banking regulatory Basel II is the sol¬ 
vency and liabilities dependence among the financial institution branches, or even 
the domino effect in the markets that could be generated by dependence among 
filial products. Thus, the solvability of each individual branch may strongly be af¬ 
fected, not only by its activities, but also by the level of dependence among all the 
branches. In consequence, it is necessary to quantify the risk, considering both the 
multivariate nature of the data and the dependence among the marginal risks. 

In Basel III, a new liquidity regulation was proposed in order to avoid the weakness 
detected in the 2007-2009 crisis; but these regulations have to be complemented 
by internal models in the institutions, in order to obtain better hedge results. These 
models have to include multivariate risk measures computable in high dimensions 
and also, to consider possible internal and external risks, even if the nature of those 
risks is strongly heterogeneous. 

In recent decades, literature devoted to extend the VaR measure to the multivariate 
setting has been published. For instance, bivariate versions have been studied in 
[Arbia (2002)], [Tibiletti (2001)], [Nappo and Spizzichino (2009)]. Also, for mul¬ 
tivariate distributions in general, some notions of VaR have been introduced (e.g. 

[Lee and Prekopa (2012), Embrechts and Puccetti (2006), Cousin and Di Bernardino (2013)]).[Embrechts and P 

linked the risk measure to the level surface defined when the distribution function 

of risk X or the survival function accumulate some a-value, which is considered as 

a quantile surface. Recently, [Cousin and Di Bernardino (2013)] introduced a new 

notion of multivariate VaR based on those level surfaces studied in [Embrechts and Puccetti (2006)]. 

They commented that considering the whole surface as a risk measure could induce 
interpretation problems. Therefore, they defined the multivariate VaR as the mean 
of the points belonging to the surface considered in [Embrechts and Puccetti (2006)] 
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and hence, the output is a point with the same dimension as the random vector of 
losses. Specifically, they define the upper-orthant Value-at-Risk ( lower-orthant 
Value-at-Risk) at a-level ((1—a)-level) as the conditional expectation of X, given 
that X stands in the o-sct of its distribution (survival) function. 

In this paper, we introduce a directional multivariate Value at Risk, based on the 
extremality level sets introduced in [Laniado et al. (2012)], which permit the con¬ 
cept of directional multivariate quantile to be defined. The extremality level sets are 
surfaces defined by following the same idea as in [Embrechts and Puccetti (2006)] 
but linked to rotations of the multivariate distribution; that is, a directional ap¬ 
proach is considered. We share with [Cousin and Di Bernardino (2013)] the idea 
that a multivariate VaR seen as a surface could bring problems in relation to its 
interpretation. Hence, we highlight the idea of considering the multivariate VaR as 
a vector-valued point that defines the vertex of an oriented orthant in the direction 
of analysis. The vertex is obtained using the mean of X to fix a reference system. 
The risk measure that we propose considers the high dimension nature of the real 
problems, and the dependence among the risks is implied in the analysis. Finally, 
we give the possibility of considering manager preferences, introducing a parame¬ 
ter of direction u. For instance, directions like the maximum variability given for 
the principal components in the portfolio, or the assets weight composition could 
be more interesting to analyze than the classic directions given for the information 
summarized in the survival or cumulative distribution functions. Besides, the direc¬ 
tional approach allows us to give bounds for the VaR related to linear combination 
of random variables, mainly when they are statistically dependent. 

We have proved properties of the directional VaR that we consider as relevant for a 
multivariate risk measure, such as consistency with respect to a particular stochas¬ 
tic order and tail subadditivity in the mean loss direction, as well as some invariance 
properties. We have compared the components of the directional multivariate VaR 
with the univariate VaR on the marginals, in order to show that the vector given by 
the VaR on the marginals provides incomplete information about the joint risk. 

We have also obtained closed expressions of the VaR when bivariate copulas are 
considered or when a multivariate Archimedean’s copulas governed the depen¬ 
dence among the components of the portfolio. Finally, we will present comparisons 
in terms of robustness with the alternative vector-valued multivariate VaR, intro¬ 
duced by [Cousin and Di Bernardino (2013)]. 

The paper is structured as follows. In Section 2, we introduce some preliminary 
concepts and notation necessary in order to understand the main contributions of 
the paper. In Section 3, the directional multivariate Value at Risk (Vai?“(X)) 
is introduced and we provide analytic properties, which can be viewed as exten¬ 
sions of those given in [Artzner et al. (1999)], to the multivariate setting. Section 
4 contains the comparisons between the univariate VaR over the marginals and the 
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components of the directional multivariate VaR. Section 5 is devoted to theoret¬ 
ical results and closed forms of the multivariate VaR when particular families of 
copulas are considered. In Section 6, we develop the robustness analysis. Finally, 
some conclusions arc outlined as well as some possible directions for future work. 

2 Preliminaries 

The main objective of this paper is to introduce a directional multivariate Value at 
Risk, based on the notion of directional multivariate quantile given in [Laniado et al. (2010)]. 
In order to make the paper self contained, we have devoted this section to revise 
the main concepts that are necessary to properly define the risk measure introduced 
in this paper. 

Definition 2.1. An oriented orthant in M n with vertex x in the direction u is defined 
as, 


= {z € M n : i2u(z-x) > 0}, (2.1) 

where u E B n (0) = {v E M n : ||v|| = 1} and R u is the orthogonal matrix such 
that R u u = e, with e = ^[1,..., 1]'. 

Based on the oriented orthant concept, we can define a partial data order (denoted 
by A u ) in M n as, 


x R u y, if and only if, (£“ D <£”, (2.2) 

where x,y Gl n . Or equivalently, 

x R u y, if and only if, i? u x < R u y, 
where the order on the right side is component-wise. 

Throughout the paper we will use the following notation related to subsets in M n . 
Given b E R r \ c E M, and A C M'\ the sets b + A and cA are defined as, 

b + 2l:={b + a:aE 2l}, cA := {ca : a E A }. (2.3) 

We recall some results on oriented orthants that will be useful in the main sections 
of the paper. 

Lemma 2.2. Given a direction u and a vertex x, then 

£x = -e=x- (2.4) 
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The proof is given in the Appendix. 

Lemma 2.3. Given c > 0 and b € M", then 

C+b = c£x + b. (2.5) 

Proof. The proof is straightforward using the definitions given in (2.3). □ 

We also recall some definitions of useful stochastic orders; see [Shaked and Shanthikumar (2007)], 
for more details. 

Definition 2.4. Given two random vectors X and Y, X is said to be smaller than 

Y in: 

(i) usual stochastic order (denoted by X < s t Y) //E[0(X)] < E[<^(Y)], for any 
increasing function <f(-) with finite expectations. 

(ii) upper orthant order (denoted by X < uo Y )ifFx(xi,...,x n ) < Fy(x \,..., x n ), 
for all x, where Fy. Fy denote the survival function ofY and Y, respectively. 

(iii) lower orthant order (denoted by X <j 0 Y) ifFy(x i,..., x n ) > Fy{xi,..., x n ), 
for all x, where I fi. Fy denote the distribution function of~X. and Y, respectively. 

It is easy to verify that both orders, the upper orthant and the lower orthant, are 
implied by the usual stochastic order. The following stochastic order defined in 
[Laniado et al. (2012)] will be a key tool in providing some properties of the mul¬ 
tivariate VaR that we will define in the next Section. 

Definition 2.5. Let X and Y be two random vectors in M n , X is said smaller than 

Y in the extremality order in the direction u (denoted by X <£ u Y)if 

E [i? u (X — z) > 0] < E [i? u (Y — z) > 0] , for all z in R n . 

It is easy to show that X <£ u Y i? u X < uo R U Y. Moreover, if X <£ u Y then 
E[X] A u E[Y], as it is proven in [[Laniado et al. (2012)], Property 3.4]. Since the 
multivariate VaR is based on the definition of a quantile, we also need to introduce 
the directional multivariate quantile given in [Laniado et al. (2010)]. 

Definition 2.6. Let X be a random vector with associated probability distribution 
function E. Then the directional multivariate quantile at level a, in direction u is 
defined as 

Qx(«, u) := 9{x € R n : P(<£) < a}, (2.6) 

with 0 < a < 1. 

From now on, we will focus on an absolutely-continuous random vector X (with 
respect to the Lebesgue measure v on M n ) with increasing marginal distribution 
functions and such that E[2Q] < oo, for i = 1, ...,n. These conditions will be 
called regularity conditions. 
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3 Directional Multivariate Value at Risk 


In the univariate setting, the relationship between the quantiles related to the loss 
distribution and the VaR is obvious. In this Section, we propose a definition of 
multivariate VaR for a portfolio of n-dependent risks, linked with the directional 
multivariate quantile defined in (2.6). Besides, the output is a point in M n ; that 
is, a vector of the same dimension as the considered portfolio of risks. Specifi¬ 
cally, as in the univariate case, this point defines the vertex of an oriented orthant 
that accumulates a probability a, but in the direction that the investor or the risk 
management considers more convenient. 

Definition 3.1. Let ~Kbe a random vector satisfying the regularity conditions and 
0 < a < 1. Then the directional multivariate Value at Risk of~K. in direction u at 
probability level a is given by 

VaR"(X) = (Qx(«, u) f|{Au + E[X]}) , (3.1) 

where A € M. 

We must highlight that given a direction u, the I/ai2“(X) is the intersection be¬ 
tween the directional quantile at level a, and the line defined by both the direction 
u and the mean of X. We want to point out that the centrality tool chosen, the 
mean, will represent a central reference point for the random vector space, i.e., for 
the support of the associated probability distribution. As we will demostrate, the 
choice of the mean in the definition of (3.1) allows us to derive desirable and in¬ 
terpretable analytic properties related to the risk measure. However, other options 
as central reference point are possible; for example the median seen as the deepest 
point associated with a multivariate depth measure, which may provide a more ro¬ 
bust risk measure (e.g. [Zuo and Serfling (2000), Cascos et al. (201 1)]). 

To illustrate this concept, you can see in Figure 1 some examples of the risk mea¬ 
sure defined in (3.1), for three different bivariate distributions in the direction —e 
with a = 0.7. This direction makes reference to the analysis of the distribution 
function of X. Figure 2 presents examples with the same bivariate distributions, 
but in the direction e and for a = 0.3; that is, taking into account the information 
given by the survival function of X. We call these two directions classical direc¬ 
tions, but the aim of this work is to show that it could be interesting to consider 
other directions in the analysis of risk. 

Observe that in the figures, the line in direction u crossing the mean in green is 
displayed while the quantile curve is displayed in red. The VaR that we propose 
is just the intersection between the line and the quantile curve. On the other hand, 


6 


the points in blue are the points "below" the level of risk a in the corresponding 
direction; meanwhile the black points are those "exceeding" the level risk. Observe 
Figure 1, if you take any point on the blue region as a vertex of an oriented orthant 
in direction —e, then the probability of that orthant will be greater than a. It will be 
equal to a or smaller than a if the point is taken from the red line or black region, 
respectively. The same conclusion can be drawn from Figure 2 but in direction e. 





(A) Bivariate Uniform 


(B) Bivariate Exponential (C) Bivariate Normal 
Figure 1: Ua-R^^X) 



(A) Bivariate Uniform 


(B) Bivariate Exponential (C) Bivariate Normal 
Figure 2: UaRo 3 (X) 




It is desirable that the classical univariate VaR agrees with our definition of VaR in 
the case n = 1; this fact will be seen in the following; remember that the univariate 
VaR is defined as, 


VaRi- a (X) = inf{x € R : P[X > x] < a}, (3.2) 

where 1 — a is usually considered closed to 1. Moreover, the VaR may also be 
defined in terms of the distribution function as, 

VaRi- a (X) = inf{x € R : P[X < x] > 1 — a}. (3.3) 
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As P[2f < x] = 1 — f'\X > x} in the univariate setting under regularity conditions, 
then (3.2) and (3.3) are the same. To be consistent with the univariate VaR, our 
definition of multivaritate VaR agrees with the classical definition for n = 1. That 
is, we have that in terms of VaR^(X), 

VaRl(X) = V aR\- a {X) = VaR^ljX), 

where VaR^X) is related to definition (3.2) and VaRf^ 0 (X) is related to def¬ 
inition (3.3). However, this fact does not hold in the multivariate context where 
F(x) + F(x) = 1 is not true in general, being 

F(x) =P[£~ e ] =P[X<x], (3.4) 

F(x) =P[<££] =P[X> x], (3.5) 

The remainder of this section is devoted to providing some properties of VaR^(X) 

which are similar to those properties considered in the risk literature; (see [Artzner et al. (1999), 

Burgert and Ruschendorf (2006), Cardin and Pagani (2010), Rachev et al. (2008), 

Cascos and Molchanov (2007), Cascos and Molchanov (2013)]). Specifically, we 
provide properties of the multivariate VaR^(X) in terms of the [Artzner et al. (1999)] ’s 
properties related to coherent risk measures in the univariate setting. Besides, we 
have explored other properties inherent to the multivariate response such as invari¬ 
ance under orthogonal transformations. All the proof for the following results is 
given in the Appendix. 

Property 3.2 (Non-Negative Loading). If A > 0 in (3.1 ), then 

E[X] A u VaR"(X). (3.6) 

This property reflects that the risk measure is a bound of the mean value of the 
losses, with respect to the partial order given in 2.2. Note that the hypothesis A > 0 
is necessary, especially when a is chosen to be close to 0. 

Property 3.3 (Quasi-Odd Measure). V a R''(-) holds the property: 

VaR^(-X) = —VaR~ u (X). (3.7) 

This property shows symmetry with respect to the random losses distribution. 

Property 3.4 (Positive Homogeneity and Translation Invariance). Let c € 
b G M n and Y = cX + b, then, 


VaR"( Y) = cVaR”(X) + b. 


( 3 . 8 ) 


Property 3.5 (Consistency w.r.t. extremality stochastic order). Let X and Y be 
random vectors satisfying the regularity conditions. 7jfE[Y] = cu + E[X] with 
c > 0, and X <s„ Y, then: 

VaR"(X) < u VaR"( Y). (3.9) 

Property 3.6 (Orthogonal Quasi-Invariance). Let Q be an orthogonal transforma¬ 
tion. Then, 

VaR$ u (QX) = QVaR^(X). (3.10) 

Property 3.7 (Non-Excessive Loading). Let R u be the orthogonal matrix de¬ 
scribed in (2.1). Then, 

VaR"(X) A u R' u sup{R u X(io)}. (3.11) 

This property shows that VaR^(X) is upper bounded by the supreme of the losses 
in the direction considered. Another good property which is desirable in the litera¬ 
ture for risk measures is the subadditivity. As is well-known, the classical univari¬ 
ate VaR is not a subadditivity measure. However, there are conditions that ensure 
the tail region subadditivity property (see [Artzner et al. (1999), Heyde et al. (2009), 
Danielson et al. (2013)]). In the same way, we highlighted that the VaR^(X) is 
not subadditive in general, but we will prove that this property holds under some 
conditions. A previous definition is necessary. 

Definition 3.8. A random vector X has regularity varying, with tail index (3 if there 
is a function f>(t) > 0 that is regularly varying at infinity with exponent jj and a 
non-zero measure p(-) on the Borel a—field £?([0, oo] n \{0}) such that, 

tP[(0(i)) _1 X G ■] 4 Ai(-), (3.12) 

when t —y oo (see [Jessen and Mikosh (2006), Resnick (1987)]). 

In this case, the measure has the property 

p(cB) = c~Pp(B), (3.13) 

for any c > 0 and B a Borel set. 

With this definition, we can state the tail region subadditivity property of the VaRf(-). 

Property 3.9 (Tail Region Subadditivity). Let X and Y be random vectors, with 
the same mean m. If (X, Y) is a regularly varying random vector with index 
(3 > 1 and non-degenerate tails then, the Voli'f i-) is subadditive in the tail region 
in direction u = i.e., 

VaR"(X + Y) A u VaR"(X) + VaR"( Y). (3.14) 
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Note that the Property 3.9 could be extended to random vectors with means satis¬ 
fying E[X] = cE[Y] for c > 0. As you can see, the property ensures that at least 
in the direction of the mean loss, it is useful to merge two risky activities in order 
to diversify the risk. 

4 Comparison of the univariate VaR componentwise and 
the Directional Multivariate VaR 

The aim of this section is to compare the components of VaR^(X.) with the uni¬ 
variate VaR related to each marginal distribution of X. But prior to this we need to 
remember the definition of a multivariate quasi-concave function. 

Definition 4.1. A multivariate function g : R n —)• R is a quasi-concave function if 
the upper-level set U q := {x € R n : g(x) > q} is a convex set for all q E R. Or 
equivalently, the complementary of the lower set L q := {x E R n : g(x) < q\ is a 
convex set for all q E R. 

We want to point out that both the distribution and survival functions, in general, 
satisfy Definition 4.1. This fact was proved in [Tibiletti (1995)], and therefore it is 
not a restrictive condition for the functions considered in this paper. Let us denote 
by Xj the i-th marginal of the random vector X and by [•], the i-th component 
related to a point in R n . The following result provides comparisons between the 
components of the multivariate VaR introduced in this work and the classical uni¬ 
variate VaR. 

Proposition 4.2. Consider a random vector X satisfying the regularity conditions. 
Assume that its survival function F is quasi-concave. Then, for all a E (0,1): 

VaRi- a (Xi) > [VaR^ (X)] i , for all i = l,...,n. 

Moreover, if its multivariate distribution function F is quasi-concave, then, for all 
a E (0,1), we have that 

[Yoi^ a (X)].> Yai?i_ a (Xi), for all i = l,...,n. 

The proof is given in the Appendix. As you can see, the preceding result can be 
extended considering other directions as follows. 

Corollary 4.3. Let X be a random variable satisfying the regularity conditions 
and fix a direction u. If the survival function of f? u X is a quasi-concave function, 
then, for all 0 < a < 1, 

VaRi- a ([R u X]i) > [R u VaR^ l (X .)] i , for all i = l,...,n. 
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Besides, if R u X has a quasi-concavity cumulative distribution, then 

[R U V a R\- a (X)] j > V aRi- a ([R u X]i ), for all i = 

where R u is the orthogonal transformation defined in (2.1). 

The proof is straightforward from Proposition 3.6 and Proposition 4.2. Therefore, 
by linking the previous results we have the following inequality for all pairs (u, a), 

(—u, 1 — a). 


VaR"(X) A u VaRf u Q (X). (4.1) 

This relationship allows us to define a directional upper VaR and a directional 
lower VaRin a similar way to [Embrechts and Puccetti (2006)] and [Cousin and Di Bernardino (2013)], 
but with a unified notation. Specifically, we have introduced the following defini¬ 
tions: 

The upper VaR in direction u is, 


VaR“(X) = VaR*(X), 
The lower VaR in a direction u is, 


(4.2) 


Mq U (X) = VaRf u JX). 


(4.3) 


An example of these concepts is displayed in Figure 3, where we can see in a 
bivariate normal distribution, the upper VaR in direction u = (-^=, -^=) for a level 
of risk a = 0.3, and the corresponding lower VaR in direction — u and level 
risk 1 — a. Note that we can describe in the plot types of asymptotes for the 
quantile curves, furthermore these asymptotes will be the univariate quantiles for 
each marginal of the rotated random vector /?, U X at the same a, where the rotation 
matrix R u is the same as in (2.1). These asymptotes can be seen as a generalization 
of those defined in [Belzunce et al. (2007)] for the quantile curves in the classical 
directions. 

Another practical situation where the link between the multivariate VaR and the 
univariate VaR is interesting (see e.g. [Embrechts and Puccetti (2006), Wang et al. (2013), 
Bernard et al. (2014)]), is when it is necessary to give bounds of the univariate VaR 
over a linear transformation of the marginal losses; for instance, when the trans¬ 
formation by the portfolio weights vector is considered, i.e., when the objective 
random variable is 


Z = w'X, 
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Figure 3: Lower and upper VaR^(X) with u = (-^=, -^=) and a = 0.3 for a 
bivariate Normal. 


where w is the vector of the portfolio weights. Since it is difficult to obtain the 
VaR of Z mainly when the components of the portfolio are not independent, there 
is special interest in obtaining at least a bound for VaR a (Z). Fortunately, we can 
give an upper-bound using our directional approach. 

Proposition 4.4. Let u = — ^ ie un ‘ tai J vector in direction of the portfolio 

weights. If ic G Qx(a, u), then w'x > VaR a (Z). 

The proof is given in the Appendix. 

Specifically as a consequence of Proposition 4.4, we have that 

— w 

w'VaRa 11 w| 1 (X) > VaRa (Z ). (4.4) 

This result is another justification to consider a directional approach of the multi¬ 
variate VaR, as well as its utility in financial applications. 

5 Directional multivariate VaR and copulas 

Researchers refer to copulas as "the multivariate distribution functions whose one¬ 
dimensional marginal distributions are uniform in [0,1]". For an extensive discus¬ 
sion of copulas, we refer the reader to [Nelsen (2006)]. This powerful tool allows 
the definition of scale-free measures of dependence and families of multivariate 
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distributions. Two aspects are important in multivariate distributions, the distribu¬ 
tion of the marginals and the dependence structure among them. The concept of 
copula fully describes the overall structure of dependence between the marginal 
variables and provides a global model for their stochastic behavior. The impor¬ 
tant result that links these two aspects is Sklar’s theorem that allows, in terms of a 
copula, to write the multivariate distribution function as, 


F(x i, • • • ,x n ) = C(Fi(x 1 ), • • • ,F n (x n )), (5.1) 

where F is the join distribution function, F\,.... F n its marginals distribution and 
C the copula, which according to Sklar’s theorem always exists. The copulas be¬ 
come a powerful tool to find closed expression of multivariate quantiles for special 
families of copulas. For example, in finance when the losses are modeled in per¬ 
centage terms, it is of practical importance to find closed expressions for the risk 
measures expressed in terms of the copula since the support of the losses will be 
the unitary hyper cube of dimension n. 

Hence, the objective of this section is to analyze how the VaR^(X) can be ob¬ 
tained in terms of some families of copulas. The first result shows the representa¬ 
tion of the VaR^(X) restricted to bivariate copulas. Let X be a bivariate random 
vector with marginals uniformly distributed in the interval [0,1]. In this case, the 
distribution function of X is a copula with density c(-, •). It is well known that 
E [X] = Note that assuming n = 2, a direction u = (ui,u 2 ) can be 

characterized by a angle # such that tan 0 = u 2 /ui, and then, u = (cos#, sin#). 
Following with the notation given by the angles, the VaR'* (X) must be a point on 
the line lg defined by, 


la : — 


f / \ Uii sin(0)—t(sin(0)—cos(0)) j 

{(w u w 2 ) : w 2 = -^-j 

{(w 1 ,w 2 ) : Wi € [ 0 , l],w 2 = \} , 


if cos(#) / 0, 
if cos(#) = 0. 


Therefore, given a direction #, VaR™(X) is characterized by its first component 
and the second one is obtained using (5.2). Now, the first component can be ob¬ 
tained by solving the following integral equation. 


Dg(w i) 


c(s, t)dtds = a, 


(5.3) 


where Dg(w\) is given by the intersection of the unitary square [0,1] x [0,1] and the 
oriented quadrant with direction determined by # and vertex (w i, Specifi¬ 

cally, Dq (w\ ) can be expressed in terms of the unknown w\ by using the semi-lines 
fg(u>i), lg(wi) that bound the corresponding quadrant which are defined as, 
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le(y> i) : = 


( 21 , ^ 2 ) • cos (j) — — ^ — z\ sin (^6 — — ^ = w\ ^tan(0) cos (j) — — ^ — sin ^6 — —^ ^ — — (tan(0) — 1) cos (j) — 


l%(w!) : = 


^ (zi, Z2 ) : Z2 sin (j) — — ^ + z\ cos (j) — — ^ = w\ ^tan( 0 ) sin (o — — ^ + cos if- l)) “ ^( tan ( 61 ) “ !) s in if — 0 j 

For instance, if # E (|, f), we can write the integral equation as follows: 


rw i 

J mini 


min{i|(tui) f|{zi=0},0} Jl%(w i) 


c(s, t)dtds+ 




'll( w l) 


c(s, t)dtds = a. 

(5.4) 


Figure 4 shows a case of the region Dg(w i) with 9 E (f, f) being the solution to 
(5.4), a point over the line lg. In summary, we can obtain VaR^(X) for a given 
bivariate vector with copula density c(-, •). 

Now, we will focus on the Archimedean family of copulas broadly used in the 
literature whose definition is the following: 



Figure 4: Quadrant given by 0 £ (f, f) and vertex over the line lg. 


Definition 5.1 (Archimedean Copulas). Let 0 : [0,1] —>• [0, oo) be a continuous, 
convex and strictly decreasing function with (1) = 0. Let (A 1 (•) be a pseudo¬ 
inverse function of (/>(•). Then an Archimedean copula C(v\. ■ ■ ■ , v n ) is defined 
by 
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i =i 



C(vi,--- ,v n ) = (j) 1 ((/)(v i)4- \-(j>{v n )). (5.5) 

In this case, for an n-dimensional random variable with distribution function as 
belonging to the Archimedean family of copulas with generator </>(•), VaR~ e (~K) 
is given by the vector with all components equal to 

[VaRt* a {X)]i = r 1 ( 0(1 ~ a) ) • (5-6) 

Moreover, if X has a survival copula C belonging to the Archimedean family with 
generator $(■), the equivalent Sklar’s representation gives the relation Fx(xi, ■ ■ ■ , x n ) 
C{F\{x{), ■ ■ ■ ,F n (x n )), where F is the join survival function and F\. ...,F n its 
marginal survival functions. Hence, we obtain that: 

[VaRKX.^ = 1 - Jr 1 ■ (5-7) 

Remember that if a vector X has a copula C, then the survival copula of 1 — X 
will also be C. Therefore, if X = 1 — X, then the copula of X and its survival 
copula arc the same; for example, Frank’s copula in the Archimedean family holds 
this property as well as the elliptical family of copulas. Then, in this case the closed 
expression for VaR%(X.) is the reflection point of Rai?j"_f a (X) with respect to the 
point (!,••• 

Now we will present some examples using some Archimedean copulas. Firstly, we 
are going to use Frank’s subclass to present an example of Vai?“(X) for any direc¬ 
tion u in the bivariate case. Later we will present some comparisons between the 
lower orthant VaR = VaR a QC) and the upper orthant VaR = VaR a (X.) devel¬ 
oped by [Cousin and Di Bernardino (2013)] with the VaR^(X.) but considering 
a n-dimensional copula belonging to Clayton’s subclass. Let’s define these two 
subclasses. 

(i) Frank Copula: The generated function of this copula is 

/ p-/3r _ i \ i 

= -In l _ i j and c/)- 1 ( s ) = --ln( 1 - (1 - e _/3 )e _s ), 

(5-8) 


Cp(vi,v 2 ) 
cp(v i,v 2 ) 




(e~P vl - l)(e~P V2 
e~P - 1 



/3( 1 - e -P)e~P( vl+V2 '> 

((e~P V1 - l)(e - ^ 2 - 1) - - l)) 2 ’ 


(5-9) 

(5.10) 
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where (3 € M\{0}. 


(ii) Clayton Copula: This family is generated by 


<As( r ) = ^(r p - 1) and 

</>/(«) = (1 + /3s) 1//3 , 

(5.11) 

Cp(vi,V 2 ) = max ji 

+ l) 1//3 ,0}, 

(5.12) 

where /3 £ [— 1,0) U (0. + oo]. 




In Figure 5 we have drawn the first component of the directional Vai2“(X) for a 
bivariate random vector, with density given by the Frank copula density. The left 
plot is related to u = —e and the right plot is related to u = — -^=(1,2). Both plots 
present the changes as 0 < a < 1 for different values of the dependence parameter 

13. 




a) Direction —e = —^[1,1]' b) Direction u = —§] 7 

Figure 5: Behavior for the first component in VaR^(X) varying a. 


We can see in Figure 5 the dependence of VaR^(X) with respect to 3. Note that 
as P —>= ±oo and the direction is ±e, we will get the extreme cases known as 
comonotonic and counter-mono tonic, respectively. In the left plot, it can be seen 
that the comonotonic case matches with the vector composed of the univariate 
VaR on the marginals, which in this case is given by the vector \VaR~ e {X)]i. In 
addition, it is well known that rotations over random vectors do not preserve the 
dependence structure in the rotated distribution; furthermore, this fact is captured 
in the right plot where the change of direction shows the rotations of the measure 
in each dependence parameter considered. 
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Let X be a random vector with distribution function belonging to the Clayton 
copula subclass. Hence 1 — X is a random vector with Clayton survival cop¬ 
ula. We have presented the comparison of the first component of VaR~ e (X i) 
with VaR^ (X) = E [X|F(x) = a] and VciR%{ 1 — X) with VaR a ( 1 — X) = 
E [X|F(x) = 1 — a], the correspondent lower ortliant VaR and upper orthant VaR 
developed by [Cousin and Di Bernardino (2013)]. 

Table 1 contains the explicit expressions of VaRg ^X) and VaR a ( 1 — X) in di¬ 
mension 2, and the generalized expressions for our proposal in terms of a, 3 in 
any dimension. Figure 6 shows the graphical comparison for n = 2; the left 
plot presents the results for VaR~ e (X ) in solid line and VaR a (X ) in dashed 
line, while the right plot presents the results for VaR° (1 — X) in solid line and 
VaR a ( 1 — X) in dashed line. 



Directional VaR^(-) 

[Cousin and Di Bernardino (20 13)] ’s VaR 

X 


(l+a-ey? 

K n ) 

P aP—a 

/3-1 a P~ l 

1-X 

1 - 

(l+(l-a)-P\~ij 

1 3 (l-«) /3 -(l -a) 

{ - ) 

0-1 (1—— 1 


Table 1: Clayton’s Copula Case 

The results in Figure 6 also shows us that in the case of random vectors with Clay¬ 
ton copula class, VaR~ e {X) increases with respect to the parameter a and de¬ 
creases in the parameter f3. On the other side, VaR%( 1 — X) is an increasing 
function of the parameter a, but also an increasing function of the dependence 
parameter (3. These features for this class of copulas were commented on and 
proved by [Cousin and Di Bernardino (2013)] and for our risk measure can be eas¬ 
ily proved following the same scheme. In addition, we need to highlight that for 
each fixed pair (a, (3), the following relationships hold, 


VaRj X) < VaR- e (X) and VaR %(1 - X) < VaR a {l - X), (5.13) 


where the inequalities are componentwise. Hence, we can say that our measure¬ 
ment is more conservative in the upper case and we are more optimistic in the 
lower case. This can be taken into consideration by the manager according to 



The previous section presents the analytic results for random vectors with [0,1]- 
uniform marginals distributions. However, in practical situations, it is necessary 
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a) Lower Case b) Upper Case 

Figure 6: Comparison for Clayton’s family of copulas. 


to obtain Uai?“(X) for any random vector X. In this case, we use the following 
computational approach summarized in the following pseudo-algorithm: 


Input: u, a, h and the multivariate sample X m . 
for i = 1 to m 

= [£;,]. 

If \P t — a | < h 

x * e 2x m («> u )> 

end 

for x j € Q Xm (a, u) 
dj = dist(xj, {/r Xm + Au}), 
end 
end 

VaR%(X m ) = {x fc | d k = minjdj}}, 


, x m } is the sample of the random vector X, /r, x the 


a,u : = 


{xj: 


— a 


< h 


| the sample quantile 


where X m := {xi, 
sample mean, Q Xm ( 

curve with a slack h and P Xm [■] is the empirical probability distribution of X, 

Using this procedure, we are able to deal with high dimension random vectors. We 
are aware that this procedures can be improved using more sophisticated tools of 
the non-parametric statistics, but they are material for another project. 

On the other hand, it is well known that in risk theory, it is desirable that a measure 
be robust, (see [Artzner et al. (1999), Burgert and Ruschendorf (2006), Cardin and Pagani (2010), 
Rachev et al. (2008)]). But in general, most of the measures are sensitive to atypi¬ 
cal observations. In this section, we present a simulation study in order to describe 
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the sensitivity of our proposal, using the [Cousin and Di Bernardino (201 3)] ’s mea¬ 
sure as a benchmark. 

The contamination model that we will use in the simulations is the following: 

x . d [ x i with probability p = 1 - w, ^ 

1 X 2 with probability p = uj. 

where x i = JVi(mi, Hi), x 2 = A^/Xi + A M , Hi + As) and 0 < uj < 1. The 
parameters of x i are, 

Ml = [50,50]', Ei = (^ “)• 

x i remains fixed in the analysis, but the parameters of the normal distribution of 
X 2 are changed to different steps to generate outliers. As a measure to quantify the 
effect of the outliers, we define, 

|| Measure(X. u ) — Measure(X °)\\2 
\\Measure(X9)\\ ’ 

where Measure(XP) is the risk measure evaluated in x °, with uj = 0 and Measure(X. ul ) 
is a risk measure evaluated with the sample with a level of contamination uj%. 


Scenarios 

Parameters of x 2 distribution 

Variance Analysis 

rfi i Hi + 

4.5 0 

0 6.5 


Covariance Matrix Analysis 

Hi + 

'4.5 0.2' 
0.3 6.5 


Mean Analysis 

Mi + A m , Hi 

Join Analysis 

. 4.5 0.2' 

Mi + A m , Hi + Q 3 6 5 


Table 2: Simulation Stages and Parameters 

We have considered the scenarios for x 2 , described in Table 2. The procedure is 
the following: firstly, we have generated a non-contaminated sample x w , uj = 0 
with 5000 observations and we calculate both VclRq i( x ) and Vaf?o.i( x )- 
Secondly, we have used the contamination model (6.1) taking values for uj from 
1% to 10%. Then, we generated for each uj, 5000 samples of xl with an expected 
value of outliers uj%. We have evaluated the risk measure as well as the percentage 
of variation for each level of contamination, performing this procedure 100 times 
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Figure 7: Percentage of variation of the measures varying the variances 


and we have reported the average of PV U in the following plots. The first scenario 
suggests outliers given by changes on the variance of the marginals, which are diffi¬ 
cult to detect in practice. We can see in Figure 7 that the behavior of VclRq x (X) is 
better than that corresponding to upper-VaR in [Cousin and Di Bernardino (2013)] 
for any level of contamination. "Better", in this context, means that is smaller. 
The second scenario considers changes in all the components of the covariance ma- 


♦* v 

v 






/*' 


Figure 8: Percentage of variation of the measures varying the covariance matrix 


trix. The results are reported in Figure 8 that shows again the better behaviour of 
VaR o i(X) with respect to robustness. The last scenarios consist of changes in the 
mean. Firstly, we affected the first component of the mean and then we affected 
the second one and finally both of them simultaneously. 

Figure 9 summarizes the results. As we can see, VaRg , (X) shows robustness 
under the presence of outliers of high dimension, but an extra-sensitivity under 
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outliers in a unique component. The use of the mean of the random loss as the cen¬ 
tral point in the definition of our VaR could be the cause of this lack of robustness. 



a) A^ = [0,25]' b) = [25,0]' c) A M = [25,25]' 

Figure 9: Percentage of variation of the measures varying all the parameters 


7 Conclusions 

In this paper we have defined a multivariate extension of the classical risk measure 
VaR based on a directional multivariate quantile recently introduced in the liter¬ 
ature. Specifically, we have proposed the directional multivariate Value at Risk 
(VaR™(X)) as a tool to analyze a portfolio of n heterogeneous and dependent 
risks considering external information or manager preferences. 

We have analyzed the analytic properties of VaR^(X) in the same way as the 
[Artzner et al. (1999)]’s axiomatic. We have provided some invariance properties 
as well as consistency and tail subadditivity property, which are desirable in a 
risk measure. We have shown relations between the components of the output 
of VaR^(X) with respect to the corresponding univariate VaR over the marginals. 

A link between the univariate VaR over the linear transformation using the portfo- 

_ w 

lio weights vector w, and the value of this transformation over VaR a |w| (X) is 
given. We have also presented closed expressions for I/ai2“(X) in terms of some 
families of copulas, considering particular dimensions or particular directions. 

Finally we have presented a simulation study of robustness comparing the behavior 
of VaR^(X) with respect to the risk measure proposed in [Cousin and Di Bernardino (2013)]. 
The simulations show advantages of our proposal in relation to the presence of out¬ 
liers. We have also detected in this study an open question to be taken into consid¬ 
eration in future work. The idea is to consider another central point instead of the 
mean as the center of the reference system, in order to improve the robustness of 
the risk measure, but, at the same time, keeping the good properties that have been 
proved. One option is to use a multivariate depth measure to choose the central 
point. 
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Appendix 

Proof of Lemma 2.2. Since, 

R u u = e and f?_ u (—u) = e, 
we have that R u = — R- u . Then using (7.1), we have that, 

= {z € M n : R u (z - x) > 0} 

= {z € M n : R-u {—z — (—x)) > 0} 

= -£” u 
V ~—X ' 


Proof of Property 3.3. Due to Lemma 2.2, it is easy to prove that 

Q_x(«, u) = — <2 x (cr, -u), 


and hence, 


(7.1) 


□ 

(7.2) 


Q-x(«, u) f|{Au + E[—X]} = (-Q x («, -u)) f| (-{A(-u) + E[X]}) 

= - (Qx(«, -u) P){A(— u) + E[X]}) . 

Then, 


VaR^(-X) = -VaR~ u {X). 


Proof of Property 3.4. This property is derived using Lemma 2.3. 

Proof of Property 3.5. Since X <£ u Y /i u X < uo R u Y, we get: 

L x (a, u) := {z € 1" : P x (£") <a}7{z6 R” : P Y (C") < a} := L Y (a, u) 

Besides, Yaii^(X) = <9Lx(a, u) P|{Au+E[X]} and VaR^(Y) = dL-y{a, u) P|{Au+ 
E[Y]}. Therefore, using the partial order defined in (2.2) there are three possibili¬ 
ties for s, t € M n : 
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□ □ 


(i) s t, 

(ii) s 2 < u t and t 2 < u s, 

(iii) s t. 


We can prove that the two first options are not possible for the points Vai?“(X) 
and VaR^i Y). Suppose that 

VaR"(X) y u VaR"{ Y), 


which implies that. 


"Zx 


C £ 


U 

Z Y ' 


Hence, 

f y(£z y ) > f x(C y ) > F x(CJ = «■ 

In which case we arrive at a contradiction,if we assume the regularity conditions. 
Moreover, the hypothesis E[Y] = cu + E[X], for all c > 0 and the conclusion 
E[X] E[Y] derived in [Laniado et al. (2012)] (Property 3.4.), permits us to 
reject the second possibility of ordering between the two points. Thus, the only 
option possible is, 

VaR"(X.) VaR"{Y) 


Proof of Property 3.6. First, note that: 


□ 


(A(Qu) + E[QX]} = Q{Au + E[X]}. 


Besides, we have the following relationship: 


= {z € M n : Rq u {z ~ Qx) > 0} and Rq u (Qu ) = (Rq u Q )u = e. 

Then R u = Rq u Q, which implies that = Ot”. and Pqx(£qx) = F x(^x)- 

Then, we get 


Qqx.(o',Qu) = QQx(a, u), (7.3) 

which proves the result. □ 

Proof of Property 3.7. Property 3.6 implies that, 

R u VaR"(X) = VaR e a (R u X), 
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where e = ^[1,1]'. Then, 

R u VaR"{X) < sup{i? u X(w)}, 


and the proof is complete. □ 

Proof of Property 3.9. It is easy to see that the equality in the mean implies that 
the vectors E[i? u X], E[i? u Y] and E[7? U (X+Y)] lie on the same line with direction 
vector e. Then, we can write: 


Cyaflg^x) = n[VaR e a (R u X.) hC 

= n[Yai^(.R u X)] 1 [w,oo) r \ 

where w is the vector whose components are - and [ ]i denotes the first component 
of the vector. 

For a > 0 small, - —>• oo, and then, 



a 


(RuX, RuY) € 4> 



—> y(B). 


On the other hand, we have the fact that the Borel set (^)) 1 £y aRe (R X j x 

(0, oo) n satisfies the following relation: 


-I 

a 


(R u x,R u Y) e U - U - (^afla(B u x) x (0,oo)") 


Or equivalently, 


Hence using (7.4), we have: 


[Ya^(i? u X)]! ~ ( y 
In the same way, 


i.oo) X (0,«.)“})% (I ]n. 


[lV,//"[i(„Y|]i ~ U (0,oo)"x 


, OO 


n 


0|-]n. 


(7.5) 


(7.6) 
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Now, in the case of the random variable i? u (X + Y), we have; 

^(iUX+Y)) = {(x,y) e (0,oo) 2n : (x + y) > VaR%(R u (X + Y))} 

= n[VaR*(R u (X + Y))]i • {(x,y) € (0,oo) 2n : (x + y) > w} 

= n[Ya^(i? u (X + Y))]!-C 

(7.7) 

where the inequalities in the expression are componentwise. As a consequence we 
get, 


1 . 


M | ($ {( x >y) € (0, oo) 2n : (x + y) > VaR%(R u (X. + Y))} 

Then using the last equality in (7.7), we finally get, 

[Ya«®(7? u (X + Y))]i ~ (m{{( x ,y) G (0,oo) 2n : (x + y) > w}})?</> (-) n. 


(7.8) 

Since in M n all the norms are equivalent, i.e., for two norms || • || and || • ||*, 
there are positive constants ci, C 2 such that ci|| • || < || • ||* < C 2 11 ■ ||. Then, 
whatever norm is taken, we use the transformation [[Resnick (1987)], pg. 267.+], 
x — ^ (| |x||, ||x| | _1 x) and rewrite /r(-) in terms of a new measure ?/(•) in V := {z E 
[0, oo] 2n \{0} : | |z| | = 1} as due to the property of the measure in (3.13). 

The relationship satisfying both measures for a Borel set A in V, it is given by, 


fj,(A) = / / l(r(u, v) € A) fir ^ 1+ ^drrj(du, dv). 


(7.9) 


v J o 


Then the measure of the Borel sets in (7.5), (7.6) and (7.8) can be expressed using 
II • ||i as: 


x (0,oo) n ^ 
M ^(0,oo) n x ^,°oj ^ 
/r({(x,y) € (0,oo) 2n : (x + y) > w}) 



lv 


iv 


lv 


^2 u i ) r]{du,dv), 

(7.10) 

. i / 

*52 v i) r]{du,dv), 

(7.11) 

, L / 

+ ri{du,dv ), 

, L / 

(7.12) 
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Now using the Mikowski inequality we obtain: 



(7.13) 


Hence combining (7.5), (7.6), (7.8) and (7.13), we have the result 

[Hai?®(i? u (X + Y))]i < [VaRZiRaX)^ + [VaR e a (R u Y)] 1; 

or equivalently, from Proposition 3.6 and the partial order defined in (2.2), we have 

for u = tt^ that: 

ll m ll 


VaR"(X + Y) A u VaR"(X) + VaR"( Y). (7.14) 

□ 

Proof of Proposition 4.4. By Definition 2.6, if x £ <2x(a, u), we have P[i? u (X— 
x) > 0] = a. Therefore, 

P[l , i? u (X — x) > 0] > a where 1 = [1, • • • ,1]'. (7.15) 

Since i? u u = e, we obtain, 


P[1 / 7? U (X - x) > 0] = F[^(R u u)'R U {X. - x) > 0] 

= p[V« ( x - x )>°] 

= P[w'X < w'x] = P [Z < w'x] 
Thus, (7.15) and (3.2) implies w'x > VaR a (Z). 


□ 

Proof of Proposition 4.2. The proof follows the same outline as that of [[Cousin and Di Bernardino (2013)], 
Proposition 2.4.]. Note that in direction u = e, 

^ = {zer:z> x}. 


Then we can write, 


L a = {xeK": P(£x) < «} 

= {xel n : P(X > x) < a} 

And we can assume the convexity of the [L a ] c by the quasi-concavity of the sur¬ 
vival function F, where [-] c denotes the complementary set. Now, as Qx(&, e) = 
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dL a = d[L a ] c , Vai2® (X) belongs to the set d\L a ] c . Moreover, from the definition 
of survival function we have that, 

F( oo, • • • , Xi, ■ ■ ■ , oo) > F(x) = F(x i, ■ • • , Xi, ■ ■ ■ , x n ) for all x e R" and i = 1, • • • , n. 

Then each component of a vector belonging to d[L a (e)] c is superiorly bounded by 
the univariate VaR at level p = 1 — a of the corresponding marginal. As a conse¬ 
quence, each component of VaR%(X.) is superiorly bounded by the univariate VaR 
at level p = 1 — a of the corresponding marginal and hence, the first inequality 
holds. Now for the second inequality, 

£~ e = {z € M n : z < x}. 

Then, we have, 


Li_ q = {x € M n : P(£" e ) < 1 - a} 

= {x € IT : P(X < x) < 1 - a} 

But, if F is a quasi-concave function, we have that [Li_ a ] c is a convex set and 
Qx (1 — a, — e) = <)L\ - a = 0[L\- n ] c . Therefore l / a/i > | ! _ 0 (X) belongs to the 
set [Li_ q ] c . Additionally, from the definition of distribution function, it is easy to 
show that each component of an element in [L i_ Q ] c is inferiorly bounded by the 
univariate VaR at level p = 1 — a of the corresponding marginal; hence, we obtain 
the result. □ 
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